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' Abstract. We report on high spectral resolution observa- 
. tions of five damped Lyman-O! absorbers: Zabs — 3.0248 to- 
; ward Q 0347-383, ^abs = 2.1411 and 2.8112 toward Q 0528 
-250, Zabs = 2.6184 toward Q 0913+072, and Zabs = 2.5226 
■ toward Q 1213+093. Line velocity profiles and heavy ele- 
. ment abundances are discussed. 

Nitrogen is found to have abundances less than sil- 
icon in the systems toward Q 0347-383, Q 0913+072, 
and Q 1213+093. The absorber toward Q 0913+072 is the 
most metal-deficient damped system known, with [Fe/H] 
< —3.2. The simple kinematical structure of the metal ab- 
sorptions makes this system ideal to discuss the [0/Si] and 
[N/0] ratios. We find [0/H] « -2.7 and -2.7 < [Si/H] < 
-2.2. 

By combining these data with information gathered in 
the literature, we study the kinematics of the low and high 
ionization phases in a sample of 26 damped Lyman-a sys- 
tems in the redshift range 1.17—4.38. We note a strong cor- 
relation between the velocity broadenings of the Si 11AI8O8 
and FeiiA1608 lines whatever the line optical depth, im- 
plying that the physical conditions are quite homogeneous 
in the sample. Statistically this shows that large variations 
of abundance ratios and thus large variations of depletion 
into dust grains are unlikely. The velocity broadening of 
the absorption lines, AV, is correlated with the asymme- 
try of the lines for AV < 150 km s^^. The broader the 
line the more asymmetric it is, as expected in case rota- 
tion dominates the line broadening. However this correla- 
tion does not hold for larger AV suggesting that evidence 
for rotational motions is restricted to velocity broaden- 
ings AV < 150 km s"^ The systems with > 200 
km s~^ are peculiar with kinematics consistent with ran- 
dom motions. They show sub-systems as those expected 
if the objects are in the process of merging. 



Send offprint requests to: C. Ledoux (cedric@astro.u-strasbg 
.fr). 

* Partly based on observations carried out at the European 
Southern Observatory, La Silla, Chile. 



Although the sample is not large enough to draw firm 
conclusions, there is a trend for the mean velocity broad- 
ening to decrease with redshift from 80 km s~^ at z < 2.2 
to 50 km s~i at z > 2.2. 

The kinematics of the low and high ionization species 
are found to be statistically correlated, though the high- 
ionization phase has a much more disturbed kinematical 
field than the low-ionization phase. This should be taken 
into account in any model of high redshift damped Lyman- 
a systems. 

Key words: Quasars: absorption lines - Quasars: individ- 
ual: Q 0347-383, Q 0528-250, Q 0913+072, Q 1213+093 
- Galaxies: halos 



1. Introduction 



Damped Lyman-a (hereafter DLA) systems are character- 
ized by a hydrogen column density 7V(Hi) > 10^° cm^^. 
The optical depth at the Lyman limit is large enough so 
that hydrogen is neutral. The gas is either cold (T < 1000 
K) or warm (T ^ lO"' K) for the highest or lo west column 

[19921) . As a con- 



densities respectively (Petitjean et al. 
sequence of the shape of the column density distribution, 
(Pn/dNdz oc N~~^ with /3 ~ 1.5, most of the mass is in the 
systems of highest column densities. The number density 
of the latter decreases with time presumably as a conse- 



quence of star-formation (Wolfe et al. 1986| ; Lanzetta et 



al. 1995 ). Indeed, the cosmic density of neutral hydrogen 
in DLA absorbers at z ~ 3 is similar to that of stars at 



the p resent time (e.g. Wolfe e t al. 1995 , Storrie-Lombardi 
et al. |1996| ; see also Turnshek |1997| ). 

Metallicities and dust content have been de rived from 
zinc and chromium ob servation s (Mey er et al. 1989 ; Pet- 
tini et al. |l994l |l9974 Lu et al. |l996b|) . The [Zn/Cr] ratio 



can be considered as an indicator of the presence of dust 
if it is assumed, that, as in our Galaxy, zinc traces the 



2 



Cedric Ledoux et al.: On the kinematics of damped Lyman-a systems 



Table 1. Journal of observations 









Date 


Tn qItii m PTii" 




\A/fi'\7'pl PTi ffi" n 

V V CL V ClCliti Lll 

T" n <TP f A 1 


S/N 


(I) 






f4) 


(5) 






(8) 


Q 0347-383 


17.30 


3.230 


20/08/90 


CASPEC 


13000 


4250-6100 


12 (25) 








16/12/92 


EMMI-R 


33000 


6100-9900 


12 (19) 


Q 0528-250 


17.24 


2.779 


16/01/89 


CASPEC 


13000 


3650-5100 


10 (18) 








16/12/92 


EMMI-R 


33000 


5795-9900 


13 (25) 








26/02/95 


CASPEC 


20000 


3650-5095 


13 (26) 








19/11/95 


CASPEC 


35000 


3750-4875 


14 (20) 


Q 0913+072 


17.10 


2.785 


16/01/91 


CASPEC 


13000 


4150-5600 


15 (37) 








16/12/92 


EMMI-R 


33000 


5800-9900 


13 (23) 








26/02/95 


CASPEC 


20000 


3900-5095 


12 (20) 


Q 1213+093 


17.20 


2.719 


26/02/95 


CASPEC 


20000 


3900-5120 


12 (38) 



gaseous abundances whereas chromium is heavily depleted 
into dust-grains. Although these assumptions have been 
questioned by Lu ct al. (1996b) and Prochaska & Wolfe 
(1997a), counterarguments have been given by Pettini et 
al. (1997^). The typical dust-to-gas ratio determined this 
way is of the or der of 1 /30 of the Milky Way value (Pet- 
tini et al. 1997a; see also Vladilo 1998). The corresponding 
amount of dust could bias t he ob served number density of 
DLA systems (Fall & Pei |1993| ) and solve the G-dwarf 
problem (Lanzetta et al. 1995 ). Metallicities are of the or- 
der of a tenth solar with a tendency for decrea sing metal- 
licity from z ~ 2 to z ^ 3 (Pettini et al. 1997b| ; see Boisse 
et al. 1998 for lower redshift). At any redshift however, 
the scatter is large and it may be hasardous to draw pre- 
matured conclusions from the small sample available. 



also been reviewed by Vladilo (1998). However, it has been 
shown recently that DLA systems at intermediate redshift 
are associated with galaxies of very different morpholo- 



Recently, Prochaska & Wolfe (1997b) have used Keck 
spectra of 17 DLA absorbers to investigate the kinemat- 
ics of the neutral gas using unsaturated low excitation 
transitions such as SillA1808. They show that the absorp- 
tion profiles are inconsistent with models of galactic haloes 
with random motions, spherically inf ailing gas an d slow ly 
rotating hot disks. The CDM model (Kauffmann 1996 ) is 
rejected as it produces disks with rotation velocities too 
small to account for the large observed velocity broadening 
of the absorption lines. Models of thick disks {h ^ 0.3R, 
where h is the vertical scale and R the radius) with large 
rotational velocity {v ^ 225 km s ~^) c an reproduce the 
data (see also Prochaska & Wolfe 1998| ). In a subsequent 
paper however, Haehnelt et al. (1998) use hydrodynamic 
simulations in the framework of a standard CDM cos- 
mogony to demonstrate that the absorption profiles can 
be reproduced by a mixture of rotational and random mo- 
tions in merging protogalactic clumps. The typical virial 
velocity of the halos is about 100 km s~^. 

The fact that damped systems originate in thick disks 
has been questioned previously. In particular, the metallic- 
ity distribution of DLA systems is inconsistent with that of 
stars in the thick disk of our G alaxy (Pettini ct al. 1997b ; 
see however Wolfe & Prochaska 1998 ) . Arguments in favor 
of DLA systems being associated with dwarf galaxies have 



gies (Le Brun et al. 1997). This strongly suggests that the 
objects associated with high-redshift DLA absorbers are 
progenitors of present-day galaxies of all kinds. 

The discussion of whether high redshift DLA absorp- 
tion systems are produced in large, fast rotating disks or in 
building blocks of galaxies is important since it is related 
to how present-day galaxies form, either through initial 
formation of large disks and subsequent accretion of gas 
or as a result of merging of pregalactic clumps. In this 
paper, we add observations of five DLA systems to infor- 
mation gathered in the literature to further address this 
problem. We present the observations in Sect. 2, and anal- 
yse the velocity profiles and the abundances in Sect. 3 and 
4. In Section 5, we construct an homogeneous sample of 
26 high-redshift DLA systems, and discuss the kinemat- 
ics of the low and high ionization species. We draw our 
conclusions in Sect. 6. 



2. Observations 

The data were collected at the F/8 Cassegrain focus of 
the 3.6m telescope, and at the Nasmyth focus of the 3.5m 
NTT telescope at the La Silla observatory, ESQ Chile. 
Blue spectra were obtained with the ESQ echelle spec- 
trograph (CASPEC). A 300 line mm^^ cross disperser 
was used in combination with a 31.6 line mm~^ echelle 
grating. For Q 0347-387 and Q 0913+072, the detector 
was a RCA CCD with 640x1024 pxl of 15/im square and 
a read-out noise of 25 electrons rms. The resolution was 
R = 13000. Additional data were obtained on Q 0528-250, 
Q 0913+072 and Q 1213+093 with higher spectral reso- 
lution {R = 20000; R = 35000 for Q 0528-250 as well). 
In this Tektronix CCD with 568x512 pxl of 27^m 

square and a read-out noise of 10 electrons was used. Three 
spectra were obtained through the red arm of the ESQ 
multi-mode instrument (EMMI) in the echelle spectro- 
scopic mode. A standard grism was used as cross disperser 
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Fig. 1. Velocity profiles in the normalized spectra of 
Q 0347-383: least square fit of low ionization lines of the 
z = 3.0248 DLA system. Components are indicated by a 
vertical bar. The curve below the spectrum is the noise 
rms in the adjacent continuum. 



in combination with the echelle grating #10 and a slit 
width of 2 arcsec, resulting in a resolution of i? = 33000. 
The detector was a Thomson THX31156 with 1024x1024 
pxl of 19//m square and a read-out noise of 5 electrons 
rms. For each exposure on the object, flat-field images 
and a wavelength comparison Thorium- Argon spectrum 
were recorded. The accuracy in the wavelength calibra- 
tion measured on the calibrated Thorium-Argon spectra 
is about a tenth the resolution quoted above. 

The observation log is given in Table 0. The magni- 
tudes and quasar redshifts are taken from Junkkarinen et 



al. (1991). The mean signal-to-noise ratio over the wave- 
length range considered has been computed from the pho- 
ton statistics after subtraction of the spectral lines. It is 
indicated in the last column together with its maximum 
value (inside brackets). The data were reduced using the 
ECHELLE package implemented within MIDAS, the im- 
age processing system developped at ESO. The cosmic- 
ray events have been removed in the regions between ob- 
ject spectra before extraction of the object. The exposures 
were co-added to increase the signal-to-noise ratio. Dur- 
ing this merging procedure, the cosmic-ray events affect- 
ing the object pixels were recognized and eliminated. The 
background sky spectrum was difficult to extract sepa- 
rately due to the small spacing between the orders in the 
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Fig. 2. Same as Fig. ^ for the two nitrogen line triplets 
and some other low and high ionization species in the 
z = 3.0248 DLA system toward Q 0347-383. 



blue. Instead, we have carefully fitted the zero level to the 
bottom of the numerous saturated lines in the Lya forest. 
The uncertainty on the determination is estimated to be 
5%. 

We have identified all the absorption features with 
equivalent widths larger than 5xFWHMxcr where <j is the 
noise rms in the adjacent continuum. Metal line profiles 
were fitted consistently assuming a pure turbulent broad- 
ening. Ionic column densities were derived with a least- 
squares technique using Voigt profiles convolved with the 
instrumental spread function, available in MIDAS (FIT- 
LYMAN program; Fontana & Ballester |1995| ). The oscil- 
lator strengths are taken from the compilation of Mor- 
ton ( 1991[ ) and the updated values of Savage & Sembach 
( 1996| ). The column densities and the metallicities mea- 
sured in the five DLA systems described in the next sec- 
tion are summarized in Table 2. We give the total col- 
umn densities obtained by summation over all the sub- 
components derived from the fit. Metallicities were esti- 
mated assuming that neutral and singly ionized species 
were associated with the neutral phase from which the 
H I column density originates. Upper limits of detection 
are computed under the optically thin case approxima- 
tion. When the lines are saturated, we indicate a lower 
limit obtained by assuming a mean turbulent broadening 
6 = 18 km s~^. Turbulent 6 values are given in the text 
for lines suspected to lie at the low column density end 
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Fig. 3. Same as Fig. |i| for low ionization lines in the 
z = 2.1411 DLA system toward Q 0528-250. 



of the logarithmic part of the curve of growth; otherwise 
stated, the line is optically thin and the resulting column 
density depends weakly on the Doppler parameter. 



3. Comments on individual systems 

3.1. Q 03^1-383, Zabs = 3.02^8 

The damped nature of this Lya absorber has been first 



recognized by Williger et al. (1989) with a hydrogen col- 
umn density log iV(Hi) = 20.8±0.1 cm"2. Pettini et al. 
( |1994D give log iV(Hi) = 20.7±0.1. The strongest part 
of the absorption produced by most of the low-ionization 
species is spread over less than 60 km s~^ consistently cen- 
tered at Zabs = 3.02483 (see Fig. Il|). As noted by Junkkari- 



nen et al. (1991), the existence of a strong Civ doublet 
at z = 2.385 2 acc ounting for the AA5241.0, 5249.7 fea- 
tures (Steidel 199C| ) is doubtful. As shown in Fig. [T], these 
features can be identified with O iA1302 and Si iiA1304 re- 
spectively, arising from the DLA system. Moreover neither 
C11AI334 nor SiiiA1260 lines at z - 2.3852 are detected 
down to Wobs = 0.14 A (5(7). There is also no Lya associ- 
ated absorption (Pierre et al. 1990; Wilhger et al. 1989 ). 

The ZniiA2026 and CriiA2056 absorptions associated 
with the DLA system have equivalent widths smaller than 
our detection limit, Wobs ^ 0.10 and 0.12 A r espect ively 
at the 5(7 confidence level (see also Pettini et al. 1994 ) . We 
marginaly detect C ivA1548 with Wobs = 0.24 A, while the 
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Fig. 4. Same as Fig. ^ for the sulfur line triplet and some 
other low and high ionization species in the z — 2.1411 
DLA system toward Q 0528-250. 



Si ivAA1393,1402 doublet is stronger though moderately 
weak. NvA1242 (Nv A1238 is blended) and AliiiA1854 
are undetected at the 5fT upper limits Wobs = 0.10 and 
0.08 A respectively. We do confirm the presence of Ni 
(Fig. 1^) but with a column density smaller than previously 
derived by Vladilo et al. ( |1997| ). By fitting the two Ni 
triplets together, we derive log iV(Ni) = 14.16±0.09 cm~^ 
assuming a Doppler parameter & = 19 km s^^. 



3.2. Q 0528-250, Zabs =2.1411 

This DLA system has been previously studied by Lu et al. 



(19961). We confirm the hydrogen column density derived 
by Morton et al. ( |1980D , log 7V(Hi) = 20.75 ±0.15 cm^^. 
The Fe 11 column density is mostly constrained by the 
FeiiA2374 line which, together with the other iron lines 
available, has been adjusted with a Doppler parameter b = 
20 km s~^ (see Table 2). The strongest part of the Fell ab- 
sorptions has a width of approximately 65 km s~^ (Fig. ^ 
whereas the OiA1302 and CiiA1334 lines are spread over 
160 km s^^. The Fell absorption lie at the red end of the 
C IV absorption. There is a strong Si iv component cen- 
tered at -1-100 km s^^ and spread over 105 km s~^ that is 
barely seen in C iv (Fig. ^ . For this component only, we 
derive column densities in excess of 14.55 cm^'^ and equal 
to 13.78±0.18 cm^^ for Si IV and Civ respectively. 
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Fig. 5. Same as Fig. |l| for metal lines in the z = 2.8112 
DLA system toward Q 0528-250. 
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Fig. 6. Same as Fig. |l| for low and high ionization species 
in the z = 2.6184 DLA system toward Q 0913+072. 



3.3. Q 0528-250, Zabs = 2.8112 

Most o f the metal lines have been observed by Lu et al. 
( 1996b ). Our spectrum goes further into the blue however 
and we show in Fig. ^the lines not observed previously at 
equivalent resolution together with a few others for com- 
parison. 

This system has drawn much attention because of the 
high Hi co lumn density, log A^^Hi) =21.35 cm"^ (M0ller 
& Warren 1993) and the detection of H2 molecules, log 
Af(H2) = 16.8 cm"^ at Zabs = 2.8112 (Leyshakov & Var- 
shalovich 1985^ S rianand & Petitjean 1998^ see also Cowie 
& Songaila |l995| ). The absorption rcdshift (zabs = 2.8112, 



as measured from the Niii lines by Meyer & York 1987 



and from the H2 lines) is higher than the QSO emission 
redshift (zcm = 2.77: Morton et al. |198C|) by ~3000 km s^^ 



which makes this system very peculiar eventhough the 
QSO emission redshift may have been underestimated. 
Three Lya emission-line objects have been detected within 
21 arcsec from the quasar (or 120 kpcat this red- 
shift and Qo = 0.1) by M0ller & Warren ( 1993 ) and con- 
firmed by Warren & M0ller (1996) to have redshifts within 



The low C1/H2 ratio cannot be easily explained by 
simple photo-ionization models if solar metallicity ratios 
are assumed (Srianand & Petitjean 19981 ). Indeed, the low 
Civ/Nv ratio and low abundances {Z/Zq ~ 0.2: Srianand 
& Petitjean |1998| ) indicate that the gas is most certainly 
not directly associated with the QSO (see Petitjean et al. 



1994). The fact that the three Lya objects are aligned 
on the same side of the QSO and have velocities -f 190, 
— 120, -1-110 km s~^ relative to the damped absorption for 
projected separations 9, 66 and 115 kpc respectively, 
argues against the assumption that all the gas is part of a 
large rotating disk (Warren & M0ller 1996| ). The system 
looks like a conglomera of individual clouds. The low and 
high ionization absorptions are spread over 440 km 
and are organized in two main sub-systems centered at 
-|-40 and -1-220 km s^^ relative to the H2 absorption, with 
a velocity broadening of 220 and 140 km s~^ respectively. 
It is interesting to note that the high and low ionization 
species have similar kinematical structures with well de- 
fined sub-components. Indeed, the weak Nv absorption 
is at the same redshift as the II2 absorption. This again 
argues for a group of individual clouds embedded in a 
somewhat isotropic ionizing field. 



200 km s~ from the redshift of the DLA system. The 
width of the Lya emission lines is very large (> 600 km 

s""'^) and continuum emission could be present (Warren & q 0913+072, z^hs = 2.6184 



M0ller 1996 ). This suggests that the gas is not predomi- 
nantly ionized by the quasar and that star-formation may 
occur inside the clouds. 



From the Lya line, we derive log A(II i) — 20.2 ±0.1 cm ^ 



consistent with the measurement by Pettini et al. (1997b) 
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Fig. 7. Same as Fig. |l| for low ionization lines in the 
z = 2.5226 DLA system toward Q 1213+093. 
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Fig. 8. Same as Fig. |l| for the nitrogen and sulfur line 
triplets together with the high ionization silicon lines of 
the z = 2.5226 DLA system toward Q 1213-^093. 



The low-ionization metal absorptions O i and Si ii show a 
single and weak component. A second component could be 
present in the red wing of C IIA1334 (see Fig. | in which 
the higher resolution data are used for Oi and Sin). 

The fit of all low-ionization lines together gives a best 
value 6 = 7 km s^^. This means that the lines are barely 
resolved. The true value of 6 could thus be smaller and the 
column densities larger than those indicated in Table 2. 
However, the non-detection of Si iiA1304 puts strong con- 
straints on h. Indeed in order for this line to be unseen in 
our data, the Si ii column density cannot be larger than 
13.58 whatever the b value is. For smaller h values, the 
fit would be inconsistent with the stronger Si ii lines. The 
lower limit for the Sill column density is 10^^'^ cm^^. If 
it is assumed that the above b value is the same for O i 
then the O i column density is well determined. Although 
the data are consistent with a single Oi, Sill component, 
we cannot completely rule out the possibility that the Si ii 
line, in reality, spans a slightly larger velocity range than 
O I which would imply a smaller b value for O i than for 
Sin. 

It is important to note that the FeiiA2382 line is not 
detected down to Wobs < 0.15 A at the 5ct confidence level. 

We do not detect high ionization lines (Sim, Siiv and 
C iv) at the velocities of the low ionization lines. Rather, 
the high ionization absorptions lie at -1-21 and -|-65 km s"^ 
relatively to the low ion absorptions. 



3.5. Q 1213+093, Zabs = 2.5226 



We confirm the damped nature of this system with a hy- 
drogen column density log A^(Hi) = 20.1±0.1 cm^^. The 
redshift derived from the SiiiA13G4 line is Zabs = 2.5226 
(the strongest part of the absorption is spread over less 
than 60 km s~ , see Fig.|). There is an additional satellite 
component detected by Sill absorption at +145 km s~^. 
This satellite system is however barely seen in OiA1302 
and not seen at all in CiiA1334. The nitrogen triplets are 
heavily blended (see Fig. ||) but the NiA1199 line, which 
leads to a low column density, log 7V(Ni) = 13.86+0.06 
cm~^, for 6 = 18 km s~^. Down to Wobs = 0.12 A (5(t), 
we do not detect the NvA1242 line associated with the 
possib le NvA1238 absorption mentioned by Sargent et al. 
( 119881 ). 



4. Metallicity 

The metal and dust contents of DLA systems have been 
estimated in the past few yea rs us ing weak lines of zin c 
and chromium (Meyer et al. 1989; Pettini et al. 1994). 



The ratio [Cr/Zn] is assumed to be an indicator of the 
depletion of refractory elements into dust as zinc is not 
depleted in the ISM of our Galaxy whereas chromium is. 

Some controversy has arisen recently about the pres- 
ence of dust in DLA systems. Lu et al. ( 1996b ) have ar- 
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Table 2. Metal abundances 



0347-383 (z = 3.02) 



; 0528-250 (2 
log N 
(cm-^) 
(4) 



= 2.14) Q 0528-250 {z = 2.81) 

[Z/H]* log JV [Z/H]* 
(cm-^) 

(5) (6) (7) 

~. 21.35±0.10'' 

>-1.96 >15.58'^ >-2.32 
>15.00'^ 

>15.60 

13.99±0.04'^ 

>-2.00 >15.94'^ >-2.28 

<15.37 

>14.61 >-1.72 <15.88'^ <-1.05 

>13.46'' >-1.77 >14.20'= >-1.63 

12.77±0.08^ >14.07'= 

15.26±0.04^ -1.04 16.05±0.10 -0.85 

>14.52'^ 

<-0.94 15.59±0.03'^ -1.03 

<14.53 <-1.24 

13.10±0.04= -1.33 13.65±0.12'^ -1.38 

12.38±0.10^ -1.90 

14.68±0.07 -1.58 15.45±0.11'^ -1.41 

13.22±0.06= -1.78 13.89±0.03'^ -1.71 

<12.28'^ <-1.12 13.09±0.07'' -0.91 



Q 0913-f072 {z = 2.61) 

log JV [Z/H]* 
(cm-^) 

(8) (9) 

20.20±0.10 

>14.08 >-2.67 

13.54±0.11 

<13.28 <-2.89 

<13.70 

<12.98 

14.28±0.14 -2.79 

<11.94 <-2.74 

<12.27 

13.34±0.24 -2.41 

<13.98 

13.29±0.11 

<13.92 <-1.55 

<12.53'' <-1.35 

<12.27 <-1.46 

<12.43 <-3.28 

<13.02 <-1.43 

<11.90'' <-0.95 



Ion 


log N 
(cm-^) 


[Z/H]* 


(1) 


(2) 


(3) 


Hi 


20.70±0.10'' 




Cn 


>14.87 


>-2.38 


Civ 


13.15±0.13 




Ni 


14.16±0.09 


-2.51 


Nil 


<14.73 




Nv 


<13.35 




Oi 


>15.16 


>-2.41 


VI 






Mgii 












Alii 


>13.13 


>-2.05 


Aim 


<12.04 




Sin 


>14.41 


>-1.84 


Sim 


<13.58 




Si IV 


13.45±0.12 




Sii 


14.46±0.04 


-1.51 


An 






Cr II 


<12.86 


<-1.52 


Mnii 






Fell 


>14.25 


>-1.96 


Niii 


<12.94 


<-2.01 


Znii 


<12.15 


<-1.20 



1213-1-093 (z = 2.52) 
log N [Z/H]* 



(10) 



(11) 



20.75±0.15 
>15.34 
>14.75t 



> 15.62 



>14.46t 
<15.08 



20.10±0.10 
>15.36 



13.86±0.06 



<13.82 
>15.59 



>14.74 
>15.88 

13.50±0.16 
<14.22 



<14.61 
<13.02 



>-1.29 
-2.21 



>-1.38 



>-0.91 



<-1.15 



<-1.00 
<-1.33 



* Metallicity relative t o the solar value from Sa vage fc Sembach 
tions); " Pettini et al. [1994 '' M0Uer & Warren |l993|; " Lu et al. 



( 19961) 



1996a 



Fit of the blu e comp onent only (see text for explana- 
Pettini et al. |l997b 



gued that nucleosynthesis alone can explain the element 
abundance ratios observed in DLA systems and that the 
presence of dust is thus questionable. They claim that the 
relative abundances are consistent with the bulk of heavy 
eleme nts bein g produced by Type II supernovae. Pettini 
et al. (|l997h| ) have shown however that, even though the 
nucleosynthesis history of the gas must play a role, the de- 
pletion levels of the refractory elements indicate the pres- 
ence of dust. The s ubseq uent level of extinction (Kulkarni 
et al. 1997; Vladilo 199^ ) is in qualitative agree ment with 
the direct extinction measurement by Pel et al. ( |l991 ). We 
discuss here in turn the three most interesting absorbers 
toward Q 0347-383, Q 0913-^072 and Q 1213+093. 

In the ISM of our Galaxy, sulfur i s not heavily depleted 
into dust grains (Sembach & Savage 1996 ). The gas-phase 
metallicity in the Zabs = 3.0248 system toward Q 0347-383 
is thus of the order of —1.5. The silicon and iron abun- 
dances are of the same order of magnitude, [Si/Fe] ~ 
and [S/Fe] ~ 0.5. This is consistent with solar relative 
abundances and little depletion into dust. Nickel seems 
to be slighly depleted as compared to iron. It is inter- 
esting to note however that, in our Galaxy, stars with 

[Fe/H ] 1, have [Si/Fe] - [S/Fe] ~ 0.4 (e .g. Fr angois 

1987| ) and [Ni/Fc] - (Gratton & Sneden |l99lD . This 
suggests that indeed the nucleosynthesis history of the 
damped systems may be different from that of our Galaxy 
and that the presence of dust at such a low level may be 
difficult to disentangle from nucleosynthesis effects. 

The most important observation in this system comes 
from N I. The nitrogen to silicon abundance ratio is smaller 
than solar ([N/Si] < —0.65). A similar abundance pat- 
tern is observed in the Zabs = 2.6184 absorber toward 



Q 0913-^072 (see below), and in the Zabs = 2.5226 ab- 
sorber toward Q 1213-1-093 (see Table 2). Such low values 
of the relative nitrogen abundance have been observed as 
well at Zabs = 2.2794, 2.309, 2.5379 and 2.8443 toward 
Q 2348-147, Q 0100-1-130, Q 23 44-H24 and Q 1946 -1-769 
respe ctivel y (Pettini et al. 1995| ; Molaro et al. 1998| ; Lip- 



1995t Lu et al. |l995| , see also Lu et al. |l99^ ). It is 



important to recall that the ionization correction factor 
for nitrogen is always close to unity for log iV(II i) > 20 



cm-^ (Viegas |1995|) 



The system toward Q 0913-1-072 is characterized by 
low metallicities and simple kinematical structure. This 
implies that absorptions from usually heavily saturated 
lines (O iA1302 and C iiA1334) can be used to derive abun- 
dances. The data are consistent with [C/H] « [O/H] w 
—2.7 although the fit of the lines gives a Doppler parame- 
ter 6 ~ 7 km s~^, which indicates that the lines are barely 
resolved in our spectrum and may be slighly saturated. 
We have argued in Sect. 3.4 that -2.7 < [Si/H] < -2.2 
which implies that [Si/0] might be slighly over-solar. The 
iron abundance is surprizingly low, [Fe/H] < —3.2, which, 
if confirmed, would be the smallest value amongst known 
DLA systems. Moreover, [Fe/Si] < —0.8. The abundance 
pattern is similar to that observed in the Zabs = 2.076 



system toward Q 2206-199 (Prochaska & Wolfe |l997a|) in 

which [Si/H] 2.2 and [Fe/Si] = -0.4. In the framework 

of chemical evolution models, and assuming no depletion 
into dust grains, these numbers can be explained if the 
system is in the early phase of a burst of star-formation. 
In such models, however, the oxygen abundan ce is ex- 
pected to be twice that of silicon (Matteucci et al. 1997). It 
would be interesting to investigate whether such an abun- 
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dance pattern is similar to that in Lya forest clouds with 
log 7V(Hi) - 15 cni-2 (Cowie et al. |1995D . 

The [N/Si] ratio is found to be less than solar. Given 
the low absolute metallicity, this system should be stud- 
ied in more detail to investigate the nucleosynthesis his- 
tory of DLA systems in the framework presented by Lu 
et al. ( 119981 ). Indeed, although the resolution of our data 
is slighly too low to definitively conclude, it seems that 
[0/Si] may be slighly over-solar. 



5. The metal line profiles of DLA systems 



Recently, Prochaska & Wolfe (1997b) have studied the 
velocity profiles of unsaturated low ionization metal tran- 
sitions (mostly SiiiA1808) of 17 DLA systems. They con- 
clude that the observed structures are consistent with a 
model of rapidly rotating, thick disk and inconsistent with 
disk models from a CDM struct ure fo rmation scenario 
(Kauffmann 1996 ). Haehnelt et al. ( |l998|) have shown how- 
ever that the data can be reconciled with such models if a 
detailed treatment of the gas dynamics is included. Here 
we want to address the problem of the kinematics with- 
out any a priori model in mind. We emphasize that large 
velocity broadenings arise most often in peculiar systems 
and that high and low-ionization species show some cor- 
relation in their kinematics. 

We have compiled a sample of 26 damped Lya sys- 
tems from the literature (see Table ^ for references). For 
each of them we have measured the velocity broadening 
of the profiles AF for the low and high ionization species 
at different optical depths. AT^J^^ is the velocity width 
of a low-excitation line measured as the velocity separa- 
tion of the wavelengths at which the residual in the nor- 
malized spectrum is r. We consider r = 0.9 and 0.5. Ex- 
cept for a few cases, there is a redshift at which the opti- 
cal depth of unsaturated low-excitation transitions (Znii, 
Fe ll) reaches a maximum. We take this redshift as the ori- 
gin of velocities. Vi (> 0) and V2 (< 0) are the largest and 
smallest velocities for which the residual in the normal- 
ized spectrum is r. We define an asymmetry parameter 
77 = + V2I / (Vi — V2) that is close to zero for a sym- 
metric profile and close to one for a one-sided profile. The 
results are presented in Table |3[ We emphasize that the 
velocity interval AT^ considered here is larger than that 
of Prochaska & Wolfe ( 1997b| ) . These authors have chosen 
to ignore 5% of the integrated optical depth at each edge 
of the profiles to correct for the internal broadening of the 
components. For high level of asymmetry, our definition 
of rj leads consequently to slightly smaller values of this 
parameter. After rescaling however, the distributions of 
the asymmetry parameter considered in both paper are 
similar (see Sect. U). In Column 12 of Table |, we give 
the measurements of d;^^ which is the velocity difference 
between the maximum of the optical depth for the low 
and high-ionization species respectively. 



5.1. The FeuXX1608,2586 and Si uX 1808 profiles 

The velocity broadening of the Fe iiAA1608,2586 (at r = 0.5 
or ~ 0.7) and SiiiA1808 (at r = 0.9 or - 0.1) line 
profiles are compared in Fig. ^ It can be seen that they 
correlate very closely. This implies that we can use the 
iron velocity interval that is easier to measure because of 
larger optical depths. The iron lines are observed in a large 
number of systems and we are able to compare velocities 



measured with the same ion. Prochaska & Wolfe (1997b) 
made used of several transitions because SiiiA1808 hue is 
not detected in all the systems. 

Since the correlation is very good, it implies that the 
physical conditions are quite homogeneous in the sample. 
First we can be confident that the criteria are well choosen 
to probe the densest regions where iron and silicon are 
both predominantly singly ionized. If it were not the case 
we would expect large variations in the A^(Si ii)/A^(Fe li) 
ratio as Si 11 persists much longer than Fe 11 when the H I 
column density decreases. Indeed the Af(Siii)/A^(Feii) ra- 
tio can be larger than ten for log A^(Hi) ~ 18 cm~^ (see 
Petitjean et al. 1992). It also suggests that large variations 



of abundance ratios and/or large depletion into dust grains 
are unlikely. These conclusions are consistent with detailed 
studies of DLA systems (Lu et al. 1996b[ Prochaska & 



Wolfe 1996, Pettini et al. 1997b, Vladilo 1998). We can in- 



fer from this correlation that, statistically, the local optical 
depth of Fe iiAA1608,2586 is approximately seven times 
larger than the local optical depth of SiiiA1808. Hence, 
using the oscillator strengths given by Savage & Sembach 
( 1996| ) , the Si 11 column density is on average three times 
larger than the Fell column density. The solar metallici- 
ties of iron and silicon are similar which implies that in 
these systems and on average, [Si/Fe] ~ 0.5 in the gas 
phase. 



5.2. The Fell velocity broadening distribution 

Fig. |l^ (upper panels) shows the histograms of the veloc- 
ity broadening measured from Fell at ^0.1 (left panel) 
and Ti/ ~ 0.7 (right panel). When weak wings are consid- 
ered (ti, ~ 0.1), the distribution is nearly uniform from 20 
to 250 km s~^. This is what is found by Prochaska & Wolfe 



(19971). However, when only the deepest part of the ab- 



sorption is considered (and we have given arguments show- 
ing that this part should be associated with the neutral 
phase) , the distribution is peaked at small velocity broad- 
enings. Most of the profiles have Av < 100 km s~^, while a 
few have Aw > 150 km s^^. It can be seen in Fig. |l^ (lower 
left panel) that the asymmetry parameter rj is correlated 
with the velocity broadening up to AV ~ 150 km s^^. 
For larger velocity broadenings however this correlation 
disappears. This approach co mplem ents the edge-leading 
test of Prochaska & Wolfe ( |l997b| ; Fig. 10) as it takes 
into account the correlation between asymmetry and ve- 
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Table 3. DLAS kinematics measurement 



* low 

(kms-^) 
(6) 



* low 

(kms^^) 
(7) 



(kms^^) 
(10) 



I high 



(11) 



Object 
(1) 



(2) 



log iV(Hi) 
(cm-2) 
(3) 



Ref. Transition 



(4) 



(5) 



Vlou 



(8) 



Transition 



(9) 



di«h 
(kms-^) 
(12) 



Q 0450-132 
Q 0449-134 
Q 0935+417 
Q 1946+769 
Q 0216+080 
Q 1331+170 
Q 2206-199 
Q 0201+365 
Q 2231-0015 
Q 2206-199 
Q 0528-250 
Q 2348-147 
Q 0216+080 
Q 0100+130 
Q 0201+365 
Q 1213+093 
Q 0913+072 
Q 0528-250 
Q 1425+606 
Q 1946+769 
Q 0347-383 
Q 2233+131 
Q 0000-263 
Q 2212-1626 
Q 2237-0608 
Q 1202-0725 



1.17 
1.27 
1.37 
1.74 
1.77 
1.78 
1.92 
1.95 
2.07 
2.08 
2.14 
2.28 
2.29 
2.31 
2.46 
2.52 
2.62 
2.81 
2.83 
2.84 
3.02 
3.15 
3.39 
3.66 
4.08 
4.38 



20.30 

20.00 
21.18 
20.65 
20.18 
20.56 
20.43 
20.70 
20.57 
20.45 
21.40 
20.38 
20.10 
20.20 
21.20 
20.30 
20.27 
20.50 
20.00 
21.41 
20.20 
20.48 
20.60 



1 
1 
2 
1 
1 
3 
4 
5 
1 
4 
1 
6 
1 
7 
5 
8 
8 
1 
1 
1 
8 
9 
1 
1 
1 

10 



Fe IIA2344 
FeiiA2586 
Fe IIA2586 
Fe IIA2344 
FeiiA2586 
FenA1608 
FenA1608 
FenA1608 
Fell A 1608 
FeiiA1608 
FenA1608 
SiiiA1304 
FenA1608 
FeiiA1608 
FeiiA1608 
SinA1304 
SinA1526 
FenA1608 
FenA1608 
FenA1608 
FeiiA1608 
FenA1608 
FenA1608 
FenA1608 
FenA1608 
FenA1608 



205 
85 
102 
52 
95 
95 
203 
245 
170 
15 
160 
40 
160 
77 
235 
217 
20 
400 
170 
14 
115 
235 
65 
80 
160 
110 



140 
55 
47 
34 
40 
75 
153 
205 
113 

65 
30 

95 
55 
215 
72 

325 
30 

43 

30 



0.71 
0.27 
0.49 
0.29 
0.00 
0.60 
0.70 
0.17 
0.77 

0.38 
0.00 
0.58 
0.64 
0.02 
0.03 

0.63 
0.00 

0.16 

0.00 



CivA1548 

CivA1548 
CivA1548 



C IVA1548 
CivA1548 

CivA1548 
CivA1548 
CivA1550 
SiivA1393 
SiivA1393 
CivA1550 
CivA1548 
CivA1548 
SiivA1393 
C IVA1548 
CivA1548 
C IVA1548 
CivA1548 
C IVA1548 



97 

375 
300 



62 
187 

500 
103 
360 
172 
97 
445 
445 
205 
125 
300 
280 
68 
202 
43 



1.10 

0.01 
0.50 



0.19 
0.44 

0.48 
0.26 
0.78 
0.49 
0.59 
0.21 
0.46 
0.02 
0.52 
0.53 
0.14 
0.32 
0.14 
1.09 



25 

130 
32 



7 

43 
5 

50 
80 
20 
25 
25 
225 
72 
40 
150 
50 
17 
35 
25 



REFE REN CES: (1) Lu et al . |l996b| ; (2) Meyer et al. |l995| ; (3) Wolfe |l995| ; (4) Pr ochask a & Wolfe ll997a|; (5 ) Prochaska 
Wolfe |1996|; (6) Pettini et al. |1995|; (7) Wolfe et al. [1994 (8) This work; (9) Lu et al. [l997|; (10) Lu et al. |1996a. 



locity broadening that is expected in fast rotating thick- 
disk models. 

Since it is the presence of systems with large velocity 
broadenings that mostly rules out th e CDM model (see 
Fig. 10 of Prochaska & Wolfe 1997b ), we discuss below 
these systems individually: 

• Q 0201+365 at z ^ 2462: If we strictly follow our 
procedure and measure the velocity broadening of Fell 
A1608 at r = 0.5, we find - 600 km s'^. This large ve- 
locity spread is due to a satellite of the main system, which 
has a column density log A^(Hi) ~ 19 cm~^ (Prochaska 



& Wolfe 1996). In this case, since the velocity separation 



is quite large, there is little doubt that we should con- 
sider that this satellite arises in a system that is not as- 
sociated with the main damped component. However this 
case clearly shows that the presence of satellites (as ex- 
pected in the case of merging processes) could bias the 
discussion on the amount of rotation needed to explain 
the profiles. By restricting our consideration to the main 
system, the velocity broadening is AV^ ~ 225 km s~^. The 
system does not show the edge-leading pattern. Actually 
the Lya line is fitted with 23 components (Prochaska & 



Wolfe |1996| ), none of which has log iV(Hi) > 20 cm^^ g^^^j 
four of which have log A^(IIi) > 19 cm^^. This system is 
definitively atypical. 

• Q 0201+365 at z = 1.955: The FeiiA1608 criterion 
gives Ay ^ 230 km s~^. Here again there is no obvious 
edge-leading pattern. Even the Lya line is poorly fitted 
by a damped component and the system resembles much 
more the LLS system ai z = 2.325, where strong Fell 
components are spread over 200 km s~^, whereas the Lya 
line d oes not show damped wings (see Prochaska & Wolfe 
1996|) . 



• Q 0528-250 at z ^ 2.811: The system is at Zabs > Zcm 
but is certainly not associated wit h the quasar (Warren & 
M0ller |l996t Srianand & Petitjean ^998]). It has AV^^;^^^^ = 
325 km s^^ and the profile can be decomposed in two main 
sub-systems, each with no edge -leading structure. We fol- 
low Prochaska & Wolfe ( 1997b ) and do not consider it as 
a typical damped system. 

2233+131 at z ^ 3.151: This system is pecu- 



liar for several reasons. Djorgovski et al. ( 1996 ) have de 
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Fig. 10. Upper panels: histograms of the velocity broadening of the Fe iiAA1608,2586 lines measured at an optical 
depth ^ 0.1 (upper left) and Ti, ^ 0.7 [upper right) for the whole sample of 26 DLA systems; the smaller binning 
indicates individual systems. Lower panels: asymmetry parameter (see text) versus velocity broadening for the Fell 
lines measured at Ti, ~ 0.7 [lower left) and for the high-ionization lines measured at ^ 0.1 [lower right). 



tected Lya emission from this system. The emission is red- 
shifted from the main low-ionization absorption by about 
(Lu et al. |1997D . Since the FeiiA1608 hne 
the resolution of the Keck 
^, we could only derive 

235 km s" 



270 km s 
does not reach = 
data and thus AVp^jj 



0.7 at 

^ km s" 



AT/0.9 



Lu et al. (1997) have argued that 
this system is a case for which rotation is apparent in the 
absorption profile. However it is clear from the OiA1302 
and SiiiA1526 lines that this system is composed of three 
well detached sub-systems spread over 300 km s~^ and 
centered at —270, —190 and —80 km s~^ from the Lya 
emission. There is no edge-leading pattern as those empha- 
sized by Prochaska & Wolfe ( 1997b ). Instead the two sub- 
systems at —270 and —190 km s~^ are equally strong and 
could both contribute to the Hi damping wings. It is in- 
teresting to note as well that the C iv absorption does not 
correlate with the low-ionization absorptions. The maxi- 
mum of the Civ absorption occurs at —120 km s^^, ex- 



actly in the gap between the low-ionization systems at 
— 190 and —80 km s^^. This pattern is a clear example 
of what could be expected from a disturbed object made 
up of interacting sub-units. It is qualitatively not different 
from the pattern seen in the 2;abs = 2.8112 system toward 
Q 0528-250. Rotational motions may be present in the 
sub-systems on velocity scales smaller than 100 km s^^. 

From these comments, it seems that one has to be 
very careful when discussing the kinematics of systems 
with large velocity broadenings. They often show the sub- 
systems that are expected if the objects are in the pro- 
cess of merging. Therefore the claim that the CDM model 
should be rejected has to be considered with caution. In- 
deed using hydro-simulations, Haehnelt et al. (1998) have 
reached the same conclusion. Even though these investiga- 
tions are not yet complete, the overall picture that seems 
to be emerging favors the idea that DLA systems are ag- 
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Fig. 11. Upper panels: velocity broadening of the Fell lines measured at an optical depth Ti, 0.7 versus redshift 
without considering the peculiar systems discussed in Sect. 5.2 (upper left) and velocity broadening of the CivA1548 
line versus that of the Fell lines measured at Tj/ ^0.1 (upper right); filled triangles indicate the peculiar systems at 
Zahs = 2.46 toward Q 0201+365, Zabs — 2.81 toward Q 0528-250, and Zabs — 3.15 toward Q 2233+131. Lower panels: 
histogram of the absolute velocity difference between the strongest low and high- ionization components (lower left; the 
smaller binning indicates individual systems) and Civ to Oi velocity broadening ratio versus redshift (lower right). 



gregates of dense knots with complex kinematics rather 
than ordered disks. This does not rule out however the 
possibility that part of the kinematics could be due to ro- 
tation. The correlation shown in Fig. |l^ (lower left panel) 
is suggestive of rotational motions in sub-systems on scales 
smaller than AV < 150 km s^^. 

Finally, although the sample is certainly too small to 
draw any firm conclusion, it is interesting to note that 
the mean velocity broadening of Fell at Ti, — 0.7 tends to 
decrease with redshift (see Fig. ^ upper left panel) if we 
do not include the peculiar systems mentioned in Sect. 5.2. 
If we split the sample in two sub-samples containing the 
same number of systems, the mean velocity broadening is 
80 and 50 km s~^ for z < 2.2 and z > 2.2 respectively. 
This behavior is expected in CDM models in which disc 



galaxies form as gas cools and forms stars at the centres 
of the dark matter haloes (e.g. Kauffmann 1996| ). 



5.3. The high- ionization phase 

The ex ampl es of absorption profiles simulated by Haehnelt 
et al. (199S) suggest that the comparison of the velocity 
profiles of high (e.g. C iv) and low (Si ii or Fe ii) ionization 
species may give clues on the nature of the objects asso- 
ciat ed wi th DLA systems. Although Fig. 6 of Haehnelt et 
al. (199S) is drawn for illustrative purpose, it is unclear to 
us whether the simulations are able or not to reproduce 
observed individual cases. Contrary to what happens for 
Fell at Ti, = 0.7, there is no correlation between ?7civ and 
AVciv at = 0.1 (see Fig. |l^; lower right panel). This 
is expected if the high-ionization phase has a much more 
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BJ CKm.s-1) 



Fig. 9. The velocity broadening of the FeiiAA1608,2586 
hnes, measured at an optical depth t^, 0.7, versus that 
of the SiiiA1808 line (measured at Tj, ~ 0.1) for the whole 
sample. 

disturbed kinematical field than the low-ionization phase 



(see also Turnshek et al. 1989| ) and indeed, ^V^^ is most 



of the time larger than AVp^^^^. However, there is a trend 
for AVgjy and AVp^^^ both measured at r = 0.1 to be 
correlated (Fig. |ll|). 

The velocity differences between the maximum of the 
optical depth in C iv and Fe ii is most often smaller than 
60 km s^^ (Fig. |l^). Close inspection of the spectra shows 
that the kinematics of high and low-ionization species are 
similar. Exceptions are the Zabs = 2.827 system toward 
Q 1425+606, in which the main Civ absorption corre- 
sponds to a high velocity sub-component, the Zabs = 3.151 
system toward Q 2233-1-131, in which the strongest Civ 
absorption is located in a gap between two low-ionization 
sub-components, and the Zabs = 1.776 system toward Q 
1331-1-170, in which the C iv profile is much more extended 
in velocity space than the low-ionization profiles. Such 
behaviors can b e rec ognized in the simulated spectra of 
Haehnelt et al. ( 1998 ) but should be investigated in more 
details. 

It is interesting to note that the AV(9jy/AVof ratio 
seems to decrease with increasing redshift (see Fig. |ll| ). 
A very large value for this ratio is found for the pecu- 
liar Zabs = 2.844 system toward Q 1946-f 769 which shows 
C IV absorptions symmetrically placed on both sides of the 
low ionization absorption. This general trend is consistent 
with kinematics to be more disturbed with decreasing red- 
shift. 



6. Conclusions 

We have presented high resolution data for five DLA sys- 
tems analyzing the velocity profiles and the column den- 
sities. The systems at Zabs = 3.0248 toward Q 0347-383 



and Zabs = 2.5226 toward Q 1213-1-093 have similar abun- 
dance patterns. Especially the nitrogen to silicon relative 
abundance is smaller than the solar value (which is true 
also for the system at ^abs = 2.6184 toward Q 0913-1-072). 
The system at Zabs = 2.6184 toward Q 0913+072 is the 
most metal-deficient DLA system known by now. We find 
[0/H] « -2.7 and [Fe/Si] < -0.8. A better limit on iron is 
needed. Indeed it would be difficult to reconcile this pat- 
tern with standard chemical evolution models. It is inter- 
esting to note that this system has a metal content similar 
to tha t of L yct clouds with log Af(Hi) ~ 15 (Songaila & 
Cowie |1996| ). 

From a sample of 26 DLA systems gathered from the 
literature, we show that the velocity widths of the Sill 
A1808 and FeiiA1608 absorptions are closely correlated. 
This suggests that the physical conditions in the neu- 
tral phase are quite homogeneous and that large relative 
metallicity variations and thus large depletions into dust 
grains are unlikely. We discuss the systems with Fe ii veloc- 
ity widths larger than 200 km s~^ and conclude that they 
are most of the time composed of several sub-systems. We 
argue, from the distribution of Fell velocity widths and 
the correlation between the asymmetry parameter and the 
velocity width, that rotation motions may be present in 
sub-systems on scales smaller than 150 km s~^. The ve- 
locity width of the strongest part of the Fe ii absorption 
decreases with increasing redshift. This suggests that the 
neutral regions get denser and exhibit faster motions with 
time. 
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